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ABSTRACT. In this article we prove an upper bound for the Lyapunov exponent j(E) and a 
two-sided bound for the integrated density of states N(E) at an arbitrary energy E > of 
random Schrodinger operators in one dimension. These Schrodinger operators are given by 
potentials of identical shape centered at every lattice site but with non-overlapping supports and 
with randomly varying coupling constants. Both types of bounds only involve scattering data for 
the single-site potential. They show in particular that both "/(E) and N(E) — i/E/tt decay at 
infinity at least like 1 / yE. As an example we consider the random Kronig-Penney model. 



1. Introduction 

In this article we will consider random Schrodinger operators H(u>) in L 2 (M) of the form 
(1) H(u)=H + V u , H = -^, V u = Y f a j (u)f(--j), 

where {ay(c*;)}j e z is a sequence of i.i.d. (independent, identically distributed) variables on a 
complete probability space (fi,^ 7 , P) having a common distribution measure k (i.e. F{ctj 6 
A} = «(A) for any Borel set A C K). In what follows we always suppose that k is supported 
on a compact interval and the single-site potential / is integrable with support in the interval 
[-1/2,1/2]. Moreover, the random variables are assumed to form a stationary, metrically transi- 
tive random field, i.e. there are measure preserving ergodic transformations {Tj}j e z such that 
a,j(Tku) = a.j_k{u}) for all uj 6 Q. The spectral properties of the operator ([[]) were studied 
in detail in [p|, fj[ 111 16J, 22]. The results are most complete for the case when / is the point 



interaction (see [H ]). 

The integrated density of states N(E) and the Lyapunov exponent j(E) are important quan- 
tities associated with operators of the form ((!]) (see e.g. [Q]). In particular, according to Ishii- 
Pastur-Kotani theorem Jl5| ] the set {E : j(E) = 0} is the essential support of the absolute 
continuous part of the spectral measure for H(u). 

The main idea of our approach is to approximate the operator ([]]) by means of the sequence 

n 

Hto(u)=H + £ <*(«)/(■ -i) 

j=-n 

with unchanged Hq, which converges to H(uS) in the strong resolvent sense. This differs from 
the usual approach where one puts the whole system in a box, which then tends to infinity (see 



e.g. [Ml]). In [12] (see also [1130) we used this approximation to invoke scattering theory for the 



study the spectral properties of the limiting operator (hi). Some other applications of scattering 
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theory to the study of spectral properties of such type Schrodinger operators in one dimension 
can be found in [11] and [|22|]. 

One of the important ingredients of our approach developed in [ |l2| ] is the Lifshitz-Krein 
spectral shift function. The spectral shift function naturally replaces the eigenvalue counting 
function usually used to construct the density of states for the operator ([[]). The celebrated 
Birman-Krein theorem (see e.g. [||]) relates the spectral shift function to scattering theory. In 
fact, up to a factor — 7r _1 it may be identified with the scattering phase for the pair (H( n \uj), 
H ), i.e. ^(E;uj) = —-k~ x 8^ (E; u) when E > 0, 

Here \T^ n \E)\ 2 and \R( n \E)\ 2 = \U n \E)\ 2 have the meaning of transmission and reflection 
coefficients, respectively, such that ^^(E)] 2 + \R^ n \E)\ 2 = 1. For E < the quantity 
w) equals minus the counting function for H^ n \u). 
In particular in JT2| ] we proved the almost sure existence of the limit 

(2) «*> = ...n 

n->oo 2n + 1 

which we called the spectral shift density. Also we proved the equality £(E) = Nq(E) — 
N(E), where N(E) and Nq(E) = 7r -1 [max(0, E)] 1 / 2 are the integrated density of states of 
the Hamiltonians H(u>) and H Q respectively. This result also extends to higher dimension in 
the continuous JT4{ ] and discrete cases. Also we showed that almost surely the Lyapunov 
exponent j(E) at energy E > is given as 

(3) ^ = 

n-s-oc 2n + 1 

where T^(E, (J) is the transmission amplitude for the pair of Hamiltonians {H^ n \oj), Hq) at 
energy E. We recall that *y(E) is defined as the upper Lyapunov exponent for the fundamental 
matrix at energy E of the Schrodinger operator H(u>). The connection between the Lyapunov 
exponent and the transmission coefficient \T^ n '(E)\ was recognized long ago [17, |18|]. A com- 
plete proof has appeared in Jl2|[. 



We note that the theory of the spectral shift function was also recently used to show that the 
integrated density of states is independent of the choice of boundary conditions [19] on the sides 
of a large box, in which the system is put. 

The conditions on the random variables a, and the single-site potential / stated above are 
slightly weaker than those in [12]. However the results of [12] which will be used below remain 
valid also in this more general case. 

The aim of the present paper is to prove global bounds for the Lyapunov exponent and the 
integrated density of states, i.e. bounds which hold for all E > and describe the correct 
asymptotic behavior in the limit E — > oo. These results are formulated as Theorems and 
^| below. To the best of our knowledge the first article to look for the asymptotic behavior of 
j(E) and N(E) in the limit E — ► oo is 0. The best known estimate for the integrated density 
of states is due to Kirsch and Martinelli JTo| , Corollary 3.1]. This bound however does not 
reproduce the correct asymptotic behavior of N(E) in the large energy limit. Another estimate, 
which is due to Pastur and Figotin (see [ |20[ Sec. V.ll.B]), is valid for an M-metrically transitive 
random field. Since our potential V LU (x) is a Z-metrically transitive field this estimate does not 
apply directly to the present situation. Our two-sided estimate leads to the bound (23) below 
which is very close to that of Pastur and Figotin. 

In what follows C will denote a finite positive generic constant varying with the context, but 
which depends only on / and k. 

We are indebted to Leonid Pastur for reading the preliminary version of this article. 
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2. The Lyapunov exponent 

We recall that the scattering matrix S(E) for a pair of Hamiltonians (H, Hq) on L 2 (R) at 
fixed energy .E>0isa2x2 unitary matrix 



(4) S(E) 



T(E) R{E) 
L{E) T{E) 



where L(E) and R{E) denote the left and right reflection amplitudes respectively. The transmis- 
sion amplitude T{E) can vanish only for E = (see [f§|6|]). To any S-matrix (Q) we associate 
the unimodular matrix 



HE) 



( i R(Ey 

T(E) T(E) 
L{E) 



\T(E) T (£) 

Let T a (E), R a (E), L a (E) be the elements of the S-matrix at energy E for the pair of operators 
(H + af, H ) and A a (E) the corresponding A-matrix. Also let A a (E) = U E ' 2 ' A a (E)U E /2 



with 



Explicitly we have 



U E 



A a (E) 



' e iVE o N 



R a (E) ~- 
TJW) T a (E) 

Lg(E) e -iVE 
J a (E) T(E) 



Consider the matrix 



A(E) = E{A a{ul) (E)^A a{ul) (E)} = J A a (E^ A a (E)d K (a) > 0, 

where for brevity we write a(u) instead of aij(u) with some j G Z. Let (3+(E) be the largest 
eigenvalue of A(E) and P-(E) the smallest. It will turn out below that (3+(E) > 1. Set 
j(E) = (log(3+(E))/2>0. 

The first main result of the present article is 

Theorem 1. Given the Hamiltonian (jl|) and the distribution nfor the coupling constant a, for 
all E > the resulting Lyapunov exponent satisfies the upper bound 

(5) j(E) < %E). 

In particular j(E) decays at least like 1 /yE at infinity. 

Proof. Let A( n \E; uj) denote the A-matrix for the pair (H^(uj), Hq), which by the factoriza- 
tion property can be represented in the form 



(6) Ato(E;u) = U E n - 1/2 J] K ]( .){E) ■ U-< 



-n-l/2 yj J fm rr -«-l/2 

j=-n 

In fact, this factorization property is a consequence of the multiplicativity property of the fun- 
damental matrix (see [12] for a proof and for references to earlier work). A short calculation 
gives 



(7) 



\T( n \E;oj)\- 2 = hi (A^iE^yA^iE^)) +1 
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With E denoting the expectation with respect to the measure P, by Jensen's inequality and (|7|) 
we therefore have the estimate 

e -2E{log |TW (E;u)\} < E | | T (n) w) |-2 j 

(8) = iE{tr(A( n )( J E;^)tA(")( J E; W ))} + i 

From the factorization property (Q) it follows that 

(— n n 
n^H^ n A ajH (i?) 

We will now make use of the fact that the ak{oj) are i.i.d. random variables. For this purpose 
define the 2x2 matrices Aj(E) > recursively by Aq = I and 

(10) Aj{E) = J A a (E) ^ Aj-i(E)A a (E)dK(a), 
such that in particular A(E) = A\ (E). Now it is easy to see that 

E{tr (AW^wlUW^a;))} 

(11) = tr(E(A^(E;^A^(E;u J )))=A 2n+1 (E). 



We now use the fact that the operator inequality < A < A' implies < tr A < tr A' and 
B^AB < B^A'B for all B. In particular we have A(E) < f3+(E) I from which we obtain the 
recursive estimates Aj(E) < (3 + (E)Aj-i(E) < ■ < P+(E) j I and hence 

(12) E {tr (a^ (E; w)t (E; w)) } < 2/3+ (£) 2n+1 . 
We remark that with the same arguments one proves the lower bound 

2f3-(E) 2n+1 < E (tr (A( n )(£;c;)tA( n )(£; 
The relation (Q), the estimate ( |I2| ) combined with (|J) and Fatou's lemma imply now 

1 \ogE{\T^(E;u;)\- 2 } 

ME) < - lim != 

y ' ~ 2n^oo 2n + l 

- 2«-» 2n+l 2 ; ' 

which proves the claim (|J). 

To establish the last claim of the theorem we recall the following well known estimates (see 
e.g. igl) 

(13) \T a (E)-l\ + \R a (E)\ <C J 



E 



valid for all large E > uniformly for all a in the (compact) support of k for fixed /. Using 
the estimate (JT3j) in @ gives the estimate j3 + (E) < 1 + C/y/E for all large £. Since 7(E) = 
(log (3 + {E))/2, this concludes the proof of the theorem. □ 

Since 7(i?) > 0, we obviously have the inequality (3+(E) > 1 for almost all E. We will give 
now a direct independent proof of this fact and simultaneously obtain an expression for f3 + {E). 
The matrix A(E) may be written in the form 



A(E) 



\b(E) a(E) y 
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with 

(14) a(E) = J [ lT J E)l2 ~ l) dK{a) 

(15) b(E) = - e ^y|^dK(a). 

This gives the two eigenvalues of A(E) in the form 

(16) (3±{E) = a(E) ±\b(E)\ 

Obviously a(E) > 1 and hence (3+(E) > 1. In fact, a{E) = 1 is possible if and only if 
R a (E) = for almost all a in the support of k. Then also b{E) = and f3+(E) = 1. 
Actually (if supp k has at least one non-isolated point) we do not believe there are nontrivial / 
and E for which this holds but in any case for such £"s the Lyapunov exponent vanishes as is 
easily verified (see also [12]), so this is a trivial confirmation of estimate (|5|) in this case. In the 
remaining case we trivially have f3+{E) > 1. 

As an example we consider the random Kronig - Penney model which is formally obtained 
from H(u) by replacing / with the Dirac 5-f unction at the origin. Then we have (correcting for 
a misprint on page 232 of Jl2|]) 

(17) T a (E) = (l + i ° 



2^fE 



(18) R a {E) = -i— — (l + i 



and our method still applies. This gives 

(19) a(E) = 1 + 

(20) b(E) -- 



2VE V 2VE 

(a 2 ) 



-i 



AE 

(a) (a 2 ) 



2y/E ±E ' 

Here for brevity by ( ) we denote the mean with respect to the probability measure P such that 

(a) =E{a(cj)} = / a<k(a), (a 2 ) = E{a(w) 2 } = / a 2 d K (a). 



In particular ( |l9| ) gives 



So also in this case j(E) decays at least like 1/yE as E — > oo and at least like l/E if the mean 
(a) of a vanishes, i.e. if on average the coupling constant is zero. 

3. The integrated density of states 

We denote by £ a (E) the spectral shift function for the pair (Hq + af, Hq). The second main 
result of this article is given by 

Theorem 2. For all E > the spectral shift density i{E) for the operator ([[]) satisfies the 
following two-side bound 

(22) E{£ a(fj) {E)} - r(E) < £(E) < E{Z a(fj) (E)} + r(E), 

where 

fjA ■ J 1 x f/ i^h^) u 

r{E) = mm < -, — E ' K K 



2'tt [l-\R a{ul) (E)\ 
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In particular E{£ a ( w ) (E)} and r(E) decays at least like 1 / y/~E at infinity. 
Remarks. 1. One can easily prove the following estimate 

E{£a (w) (£)} " 1 < < E{£ a(uj) (E)} + 1, 

which is valid for all B6l. 

2. By the monotonicity of the spectral shift function with respect to perturbation £(E) > if 
supp kCM + and £(E) < if supp kCL for almost all E > 0. 

3. For large E > by (fil 



K ' \2' vr \l - \R a(uj) (E)\ f f tt \l - \R a(u}) (E)\ j ~ ^E 

4. In Q we proved the relation £(E) = N (E) - N(E) = s/E/tt - N(E), where N (E) 
is the integrated density of states for the free operator Hq. Theorem ||] then gives the following 
two-sided bound for the integrated density of states 



(23) ^-E{£ a(u) (E)}-r(E)<N(E)<^-E{t; a{<Jj) (E)} + r(E), E > 0. 



TT K ' 7T 



There are some other upper bounds on the integrated density of states. A well-known result is a 
one-sided bound due to Kirsch and Martinelli JTo| , Corollary 3.1], 



C { r 1 / 2 
N(E) <-E (E + t] -V u 

y/V [J -1/2 



(x))+dx 



for any r\ > and all E S R. This bound however does not reproduce the correct asymptotic 
behavior of N(E) in the large energy limit. 

5. The bounds (^ and (22) are of interest in the context of the Thouless formula (see e.g. 



[20]) 



(24) 7 (E) - 70(E) = - j log \E - E'\ d^(E'), E G R, 

where 70(E) = [max(0, — E)} 1 / 2 is the Lyapunov exponent for H . The Thouless formula in 



the form (E4J) can be viewed as a subtracted dispersion relation (see e.g. [12]) 



Proof. In [ ]12| ] we proved (see Theorem 3.3 there and its proof) that for any two potentials V\ 



and V 2 with (compact) disjoint supports one has 

t(E; H + V 1 + V 2 , H ) = H + V U H Q ) + £(£; H + V 2 ,H ) + fa{E) 

with 

ti 2 (E) = log = — = . 

^ 12V ' 2ni & 1 - R X {E) L 2 (E) 

where R k {E) and L k (E) are the right and left reflection coefficients for the Schrodinger equa- 



tion with the potential V k , k = 1,2. Actually Theorem 3.3 in [JT2J] states that \£i 2 (E)\ < 1/2 



for all E > 0. Now we improve on this estimate. As in [12] we set 



L k {E) = a k (E)e td k , R k (E) = a k (E)e td * , k = 1,2 
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with < a,k(E) < 1. Moreover au(E) = 1 only when T^(E) = 0, which we recall can happen 
only if E = 0. Therefore 

1-R X (E) L 2 (E) 



log 



1-R X (E) L 2 (E) 



log 



l-a 1 (E)a 2 (E)e i ^ R)+s ^ 
l-a l {E)a 2 (E)e- i ^ R) + & 2 ) ) 



-2% arctan • 



a 1 (E)a 2 (E)sm(5[ R) +5 { 2 L) ] 



1 - ai (E)a 2 (E) cos(5[ R) + 5 { 2 L) ) ' 
By means of the inequality | arctan x\ < \x\ we immediately obtain 

1 1 ai (E)a 2 (E) 



(25) 



\£l 2 (E)\ <min 



2n 1 - ai(E)a 2 (E) 

Since < a k (E) < 1 we can replace a\(E)a 2 (E)(l — a\(E)a 2 (E))~ 1 either by a\(E)(l — 
ai(S))- 1 or by a 2 (E)(l -a 2 (E))-\ 

Now let us consider the operator (u>) for finite n. Applying the inequality (25) we obtain 

&\E;uj) -U( U )(E) -t a _ n{oj) (E)- ^- l \E;cj) 



< min 



1 1 \Ra n (u)(E)\ 



2 , 7rl-\R oln(u)) (E)\ j 
Repeating this procedure recursively we obtain 



+ min 



&HE;u)- £ £, J(W) (JE) 



3=-n 



< min 

j=-n 



1 1 |ik_„ I 1 

2'vrl-|i? Q _ nH (E)| J 
1 1 |flq»(£)l 1 

2' 7T 1-1^.(^(^)1 J ' 



From the existence of the spectral shift density (Q) by the Birkhoff ergodic theorem it follows 
that 

1 \R a (u>)(E)\ 11 



< E< 


f " 1 

mm < 







From the obvious inequality 



W J . Jl 1 K( W )(^)| 
E < mm < - , — - 



2'7rl-|i? QM (£)| 







||<mm|i 


H 







|-Ra(w)(-^)I 



the bound ( J22[ ) follows. 

For large E" we have the following asymptotics [||] uniformly in a on compact sets: 

a 



Ra(E) 

L a (E) 



2iy/E 
a 

2i\[E 



m 

I e~ 2i ^ t f(t)dt + 0(E- 1 ) 

m 



such that R a (E) = 0(1 /VE) and L a (E) = 0(1/ VE). If the single-site potential / has p 
derivatives in L X (E) then L a (E) = 0(E-^ P+1 ^ 2 ) and R a (E) = 0(E~ ( -P +1 ^ 2 ) as E -> oo 
[|]. The estimate E{£ a(a /,(£)} = 0(1/ VE) is Proposition | below. □ 

As an example we consider again the random Kronig-Penney model. The single-site spectral 
shift function is given in this case by 



£ Q (E) = — arctan 

7T 



a 



2VE 



E>0. 
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Therefore 



and thus 



Ek Q y(^)l = - / arctan 



a 



E 



< 



2VEJ 

<M> 



) d«( 



a 



2irVE 

Using the explicit expression for the reflection amplitude one can easily show that 



<H) + (^> <E 

2y/E 4E ~ 



R a (oj) {E) I 



We complete this section with an estimate on E{£ a ( w ) (-£')} m the general case. We will prove 
Proposition 3. There is a constant c > independent of E, f, and k such that for all E > 

C 



E 



{&(*)(£)} 



< 



2V£ 



E 



{l«Ml 1/2 }' 



1/2 
-1/2 



Let Z 1 / 2 ^ 1 ) denote the Birman-Solomyak class of measurable functions V for which 



1^1 




J+l/2 
j-1/2 



1/2' 



< oo. 



The claim of the proposition immediately follows from the following 

Lemma 4. Let V € Z 1 / 2 (L 1 ). There is a constant c\ independent ofV and E such that 



\Z(E;H + V,H )\ < 



(-i 



2VE 



\V\ 



«V2(Ll) 



for all E>0. 

Proof. As proved in [^] there is a constant C2 > independent of E and V such that 



where F 1 / 2 = signF | 1/2 , flo(«) = (iT - js) -1 , and 
(see e.g. [p]]). From the proof of Proposition 5.6 in [21] it follows that 

\\V 1/2 Ro(E + iO)\V\ 1/2 \\ Jl < 



j x denotes the trace class norm 



C3 IIT/II 



for all £ > 0. 



□ 
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